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Abstract 

A theoretical study of the improvement of the Electron Cyclotron Current Drive (ECCD) effi- 
ciency in regimes where most of the current is driven by Lower Hybrid (LH) waves is presented. A 
perturbation technique is employed to solve the adjoint equation and derive the response function 
including both collisional and LH effects in the limit where the former dominate. An alternative 
treatment of the problem, involving a numerical solution of the Langevin equations is proposed 
to gain insight into the current drive mechanism and confirm the obtained results. The existence 
of a cross-effect between the two waves is demonstrated and the conditions for the synergy, i.e. 
significant enhancement of the ECCD efficiency in the presence of LH power, are identified. 

PACS numbers: 52.35.Hr, 52.35.Mw, 52.55.Wq 



* Electronic address: Iremi . dumont @ cealrl 



1 



I. INTRODUCTION 



Among the desirable features of a future tokamak-based fusion reactor is its steady-state 
operation, which imphes that the toroidal current has to be totally sustained by non induc- 
tive sources^]. Moreover, for ongoing as well as for future experiments, a sharp tailoring of 
the current profiles is known to have a favorable effect on the plasma confinement. These 
prerequisites are key elements of the advanced tokamak concept, which relies on the fact 
that a large fraction of the plasma current is supplied by the bootstrap current, generally 
triggered and supplemented by injecting neutral beams of radio-frequency (RF) waves in 
the plasma P]. 

A wide class of waves can be launched in the plasma yet, for current drive purposes, the 
excitation of superthermal electrons has been the most successful method, in accordance 
with theoretical predictions jli]. Lower Hybrid Current Drive (LHCD) is a well-tested and 
efficient method 3, 0, O], based on Landau damping of the wave power. Its main drawback 
is that in the so-called multi-pass regime, the current profile remains difficult to control 
and calculations still lack full reliability. Note, however, that in very hot plasmas, the wave 
is expected to be absorbed without relying on any subtle upshift mechanism (single-pass 
absorption), which makes the LH wave a serious candidate to drive off-axis current in future 

n 

experiments I4|. On the other hand. Electron Cyclotron (EC) waves exploit the resonance 
with the electron gyro-motion and are characterized by a narrow deposition, thus providing 
a means to induce a local modification of the current profile. The efficiency, however, is 
known to be significantly lower than for LH waves j^. 

Owing to these complementary features, combined schemes, in which LH and EC waves 
are used together, constitute an appealing solution for advanced tokamak discharges. In 
particular, both waves are expected to exert a cross-effect on superthermal electrons, often 
and sometimes improperly referred to as the LH-EC synergy effect. The consequences of this 
cross-effect encompass the possibility to modify the LH superthermal tail, a better absorp- 
tion of EC waves and an irnprovement of the ECCD efficiency which may help in controlling 



ma an irnprc 



the current profile [7|, |8|, |9|, |10|, Self-consistent, dynamical calculations including waves 



kinetic and heat transport effects have stressed the advantages of such scenarios but owin; 



to the non-linear nature of these simulations, the various effects are difficult to separate 
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Experimentally, a cross-effect has been observed under certain conditions but the interpreta- 



2 



tion of these measurements are difficult and have led to contradictory conclusions, due either 



to their transient character or to very large fast particle losses 



a hard X-ray camera 
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16|. More recently, 



17l | has been employed to analyze the emission of fast electrons in the 



presence of LH and EC waves during the current fiat-top phase on the Tore Supra and 
FTU[19] tokamaks. In both experiments, a response of the high energy channels was ob- 
served, compatible with the existence of a cross-effect, but due to the relatively low applied 
EC power, no conclusive observation could be made in terms of ECCD efficiency. 

Despite these numerical and experimental results, and owing to the lack of a simple 
mechanism to explain the cross-effect of both waves on fast electrons, this subject is still 
being debated. The goal of this paper is thus to study the combined current drive process 
from a theoretical standpoint and identify the conditions for the existence of a LH-EC 
synergy. To this aim, the adjoint method, originally proposed by Antonsen and Chu|2(| and 
generalized to the RF current drive problem by FischQ], is extended to a situation in which 
two waves are simultaneously present in the plasma. Another possible method is to solve the 
Langevin equations, which track individual electron relaxation paths and deduce quantities 
of interest for the current drive problem by average over statistical realizations [2l|. These 
two complementary methods are employed to derive the response function, which leads to an 
estimate of the current drive efficiency, including the cross-effect and therefore the LH-EC 
synergy. 

This article is organized as follows. After the presentation of the kinetic aspect of the 
problem in section|nJ the adjoint method is employed to compute an approximate expression 
of the response function for a LHCD plasma in section IIIII Section IIVI is devoted to the 
study of the dynamics underlying the RF current drive process, through the derivation of 
the associated Langevin equations. The properties of the response function obtained from 
these two methods are examined in section The practical consequences of the presence 
of LH waves on the EC Current Drive are discussed in section IVIf where the synergy effect 
is demonstrated. Conclusions are drawn is section IVIII 
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II. KINETIC MODELING OF LH+EC CURRENT DRIVE 



In the absence of a static parallel electric field and including the effects of Coulomb 
collisions, LH and EC waves, the kinetic equation can be written as 

J--Cf = Di,f + D,J (1) 

In this expression, r = I'et is the time in terms of the collision period z/~^, / is the electron 
distribution function, u = p/ V rrieTe is the normalized momentum. 

C is the linearized collision operator. In this paper, we shall assume that for any perturbed 
distribution function written as / = + a), the high velocity collision operator Cf = 
C{fma) is linearized according to C{fma) = C{fm, fm(y) + C{fma, fm) + C{fm(y, fi), fi being 
the ion distribution function, li a = u\i/u refers to the cosine of the pitch angle, it can be 

n " 

expressed aslH 

Cf ^ -- (-^ + f)+ ^ A(i _ .2)^ (2) 
where Zi is the plasma ion charge. 

Dih (resp. Dec) is the quasilinear operator associated to LH (resp. EC) waves, which is 
related to the corresponding quasilinear diffusion tensor Dih (resp. Dec) and flux Sih (resp. 
Sec) by 

d = df d 
Dlhf = TT- ■ D;/j ■ — = — — • Sih (3) 
ou ou ou 



and 



d = df d 

Dec/=^-Dec-5^ = -^-Sec (4) 



The LH wave is absorbed in the plasma by Landau damping, which implies that Dih, the 
corresponding quasilinear diffusion tensor is dominated by its parallel-parallel component 
and the corresponding quasilinear operator can thus be written under the form 

Dih = 7^ ■ j^dih{u\\)— (5) 

Te is the local electron temperature, Dih is a constant determined by the wave power, so 
that Dihfi = Dih/ UfmeTe quantifies its relative intensity compared to colhsions. dih = dih{u\\) 
represents the shape of the diffusion coefficient. Here, we consider a regime in which the 
electrons undergo the effects of the wave in a region of velocity space bounded by two limits. 



4 



Mm and M||^2, determined by the propagation properties of the wave|12l|. The following shape 
is assumed 



,1 



Ai exp ( - - M||_,)VAMp) , M|| < M|| 
dlh{u\\) = { U\\^i/U\\ , Mm < Mil < M||,2 (6) 

Ar exp ( — (M|| — M||_r)^/ Am| 2) ) U\\ > M||^2 

Ar, Ai, u\\^r and My^^ are constants whose values are determined by the requirement that 
both dih and ddih/du\\ be continuous at My = m^ and u\\ = M||^2- 

For the problem under discussion here, a useful form for the distribution function is 
/ = + (p + where is the Maxwellian. /^(l + (p) is the distribution function 
modified by application of the LH power, solution of 

^-C'(/^0) = + (7) 

Upon subtracting ((Zj) from ((H), we obtain the equation for fmS(p, which writes 



dr 

or equivalent ly 



- C{U6<j)) = Dec/™(1 + <P + 6<j)) + Dih{fmS<P) (8) 



- c{u<p) - MfmH) = ■ s,, (9) 

In this expression, the right term describes the electron excitation caused by the elec- 
tron cyclotron waves whereas (C + Dih){fm5(j)) is representative of the relaxation under the 
combined effect of collisions and LH power. 

Recalling that is even for My, the normalized current associated to / can be written 

as 

] = j duf^cf) + 6(P) = Jo + ji (10) 

where jo and ji are 

jo = J duuiifm(j), ji = J duuiifmS(f) (11) 

Note that up to this point, no approximation other than the fact that Tg and ng, the elec- 
tron temperature and density, are non-varying on the time-scale of the studied problem has 
been introduced. To evaluate the driven current, it is usual to resort to a Fokker-Planck code 
and compute the distribution function solution of Eq. ^ In this work, however, an analytical 
approach of the problem has been preferred. Even though it implies certain assumptions, 



the linearization it is based on has the advantage of allowing a separate treatment of the 
two waves, which is impracticable in purely numerical calculations. This separation is the 
key to demonstrating unambiguously the possibility of a synergy effect. 

To compute jo (Eq. ITT|) . it is necessary to solve Eq. [7| Following Eq. 01 the expression 
for the LH-induced quasilinear flux is 

ou an 

If Sih is not significantly modified by the presence of EC waves, or equivalently assuming 
that the shape of the distribution function is mainly determined by the effects of collisions 
and LH wave, i.e. \6(f)\ ^ |0|, Eq. [7|can be rewritten as 

^-ciM = -^-s^. (13) 

Introducing go{u, u', r — r') the Green function associated with Eq. [T31 solution of 

^-Cgo = 6{u - u')S{r - r') (14) 
and defining the steady state response function as 

dr' J du'u\\go{u,u',T') (15) 



jo is given by 



jo = / duSih ■ ^ (16) 



Likewise, to compute ji, a similar treatment is applied to Eq. |H1 introducing gi{u, u', r — 
r') the solution of the Green problem for Eq. El It solves 

^ - (7(71 - Dingi = S{u - u')6{r - r') (17) 

The associated steady-state response function writes 

POO P 

Xi{u) = J dr' J duu\\gi{u,u',T') (18) 
and allows to evaluate ji by using 

ji = J duSec ■ ^ (19) 
This method allows to envision the current drive mechanism as a two steps process: 



1. A drive, whose features are contained in the expressions for the quasihnear fluxes Sih 
and Sec- 

2. A relaxation, which is described by the response functions Xo Xi- 

The quasilinear fluxes contain the information on the distribution function and their 
evaluation is a delicate task, generally involving a kinetic code. However, if one is merely 
interested in an estimation of the current drive efficiency, as long as the interaction is well 
localized in velocity space, the information on the direction of these fluxes is known to be 
sufficient . Note that this relaxes the assumption employed to derive Eq. since as a 
result, only the direction of Sih has to be unchanged by EC waves for the LH efficiency 
calculation to remain valid. 

The problem now reduces to evaluating the response function and in order to achieve 
this, two methods are available: the adjoint formalism and the Langevin equations. Both 
will be discussed for the case of two waves in the next sections. 



III. ADJOINT METHOD IN THE PRESENCE OF TWO WAVES 

In this section, the adjoint formalism is extended to the case of two waves, when no 
static electric field is present, which is relevant for a number of experiments. For the sake 
of concision, only the significant steps of this method, extensively discussed in Ref. UJ and 
references therein, will be recalled. 

By making use of Eq. El when steady state is attained, the current jo defined in Eq. ^1 
can be written as 

Jo = -j duxoC{U(j)) (20) 
Introducing the commutative operation for two functions <f{u) and ipi'^) 

[^,^] = I du<^(u)^(u) (21) 

and defining the adjoint of an operator D as 

b, Z)t^] = [bcp, ^] (22) 

Eq. EH can be rewritten as 

Jo = - / dnU^C^Xo (23) 



The adjoint equation is obtained by comparing Eqs. |221 and and making use of the 
property /^CV = C{fm^p)- It writes 

Cxo = (24) 

This equation simply describes the response of a colhsional plasma and associated with 
Eq. 121 it leads to the well-known Fisch-Boozer response function 

= WTzt/^ ^'^^ 

In order to compute ji, a similar treatment is applied to Eq. El introducing the associated 
response function Xi(u). By noting that the LH quasilinear diffusion coefficient is self- 
adjoint, i.e. -D|^ = Dih, the adjoint equation takes the form 

[C + Dih]Xi = -uii (26) 

This equation describes Xi, response function of a plasma in which the LH wave modifies 
the distribution function, modifying in turn the electron relaxation properties. Physically, it 
means that these electrons describe a collisional curve in velocity space, carrying elemental 
current u\\ + Djf^xi instead of u\\. A rigorously equivalent interpretation is that the particles 
carry u\\ as elemental current but describe relaxation curves influenced by the wave. 

Under this form, the adjoint equation does not appear to have an analytical solution. 
A further assumption is to consider that collisions dominate the relaxation process. In other 
words, in spite of the modification of the dynamics in the parallel direction caused by the 
LH wave, the relaxation curves remain mostly collisional. This approximation allows one to 
linearize xi, letting xi = X + with \6x\ <^ \x\y small parameter being Ah,o- 



The zero-th order expansion of (j26|) writes 

Cx = -uii (27) 

which demonstrates that x is exactly the Fisch-Boozer response function xo- 
To first order, we obtain 

C6x = -DihXo (28) 
which, upon expanding the quasilinear operators and letting Z = {Zi + l)/2, yields 
9Sx ^9 ^ d5x Dih 3 9 , . .dxo 



The associated Green equation writes 



Here, G^(u, u') is the steady-state Green function of the problem 
The Pitch-angle scattering term in Eq. 1201 suggests the expansion 

'^(/i' - = E ^^^^Pi{i^)W) (31) 

1=0 

where {Pi) are the Legendre polynomials and applying a variable separation, one obtains 

1=0 ^ ^ 

where Y is the Heaviside function. 
This leads to the solution for 5x 

*^(-^) - i^){£:^yti^QMmp) (33) 



with 



and 



,N Zl(l+l)+A 

Qi{u)^ I du'{-] Mu') (34) 







u 



3di,{u', /iO/i'(3 + fi'') + ^w'(3/i'' + 1) 



(35) 



IV. LANGEVIN EQUATIONS 

;hs, a natural and convenient 



To track the electrons trajectories on their relaxation pat 

25[. Besides providing a clear 



method consists in solving the Langevin equations 

insight in the dynamics underlying the relaxation process can be used to compute 

the response function x- Another advantage is that no specific approximation regarding the 
respective intensities of the collisions and of the wave has to be introduced. They can thus 
be used to validate the results obtained with the adjoint method (see section ITTH) . 

In this section, the discussion will focus on the evaluation of Xi but it can be easily 
transposed to Xo since the latter is a particular case of the former (with -D/^^o = 0). 

The Green function gi corresponding to Eq. ^ has been introduced in section |n] Physi- 
cally, gi{u, u', r)du is the probability of finding an electron initially at velocity space position 



u' within element du at location u after a time r. It means that the associated steady-state 
response function whose definition is given by Eq. UHl can be evaluated by computing 

the elemental current carried by each electron of a set whose initial location is u along its 
relaxation trajectory, as it undergoes the effects of Coulomb collisions and LH wave power, 
and perform an ensemble average afterwards. Rather than directly solving Eq. El it is thus 
possible to resort to a stochastic description of this relaxation process. It can be done by 
casting this equation into the form 



dgi 
dr 



ou 



(36) 



Introducing the friction vector F and diffusion tensor D and using Einstein convention 
for repeated indices, the probability current components are written as 



Si = FiQi 



d 



DijQi 



(37) 



By identification of Eqs. EElandEl the diffusion tensor can be written as D = Ticoii + D«/i 
with 

2/1 



D 



coll 



and 



dihiu) 



(Z, + l)(l-/i2)/2 



U/i(l — /i^) (1 — /i 



2\2 



(38) 



(39) 



For the force term, F = ^coii + ^ih with 



coll 



.,,3 



1 + — ]ueu + {Zi + l)ne^, 



(40) 



and are the unit vectors corresponding to the u and /i directions. 
To compute the drift caused by the wave, it is necessary to bear in mind that when the 
transformation from one coordinate system to another is not linear, as is the case when 
the LH quasilinear diffusion coefficient is transformed from {u±,u\\) to {u,fi) coordinates, 
the friction term needs to include a contribution from the diffusion in the first coordinate 
system j2l|. This yields 



D 



lh,0 



Ih 



ufid'if^{u) + (1 - ^i^)dih{u) Me„ + ud'if^{u) - 3fidihiu) (1 - fi'^)e 



(41) 
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with d'lf^ = ddih/du\\. 

The Langevin equations describing the electrons trajectories can then be written as 



^ = h(u,r)+i(u,r)-C(r) (42) 

The components of the matrix g are hnked to the diffusion tensor D through 

gi, = (d'^\, (43) 

= 1/2 

where D is obtained by diagonahzing D, taking the positive square root of the eigen- 
values, and transforming the diagonal matrix back. 

In the framework of Stratonovitch calculus, the proper deterministic force components 
are given bv|2l| 

= 1/2 pi =1/2 

/i, = F,-(D )fc,-^(D ),,■ (44) 

OUk 

The second term of the left-hand side of this equation is a correction to the noise-induced 
drift. Albeit straightforward, the computation of h and g is tedious, the resulting expressions 
cumbersome and they will not be presented here. 

£(r) is the Gaussian-distributed Langevin force, described by its stochastic properties 

(e.(r)) = 0, and (6(r)e,(r')) = 25,/(r - r') (45) 

where (■) refers to the average performed over statistical realizations. 

Practically, in the simulations presented in this paper, the collision time is split in numer- 
ous timesteps and each electron velocity evolves according to Eq. until thermalization is 
attained. It should be emphasized that the integration of a stochastic equation in the case 
of a multiplicative noise has to be carried out with care to avoid numerical artifacts liable 
to distort the result. A detailed discussion of this question can be found in Ref. |26|. Here, 
both Euler and Heun methods have been implemented and the results have been found to 
be generally indistinguishable. It is recognized that two main sources of numerical error can 
alter the solution of stochastic equations^]. Firstly, the statistical error, due to the finite 
number of realizations, which can be evaluated from standard statistical methods. Secondly, 
the error induced by the time discretization. To reduce the latter, it is usual to perform 
several simulations with different timesteps Ar, and extrapolate the result for Ar 0. For 
the values of Dihfl considered in this section, this error was found to be negligible. 
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To illustrate the combined effects of Coulomb collisions and LH waves on the electron 
relaxation, randomly chosen individual relaxation paths can be studied. The velocity-space 
configuration of the problem appears on Fig. Q where three domains (labeled 1-3) are 
distinguished, corresponding to (1) Uq < u\\^i, (2) uq > u\\^i andMyo < u\\^2, and (3) up > u\i^2- 

Such sample trajectories are shown on Figs. |21 El IH and which have been produced 
with the following parameters: = 3, u\\^2 = 5, Au\\^i = 0.5, A'U||^2 = Zi = 1 and 
Dihfi = 0.1. All these figures are divided in two, the part labelled (a) where only the 
collisions are included in the calculation, and (b) where the electrons undergo the combined 
effect of collisions and wave power. 

On Fig. 121 a trajectory obtained for an electron whose initial velocity lies in region 1 is 
shown. In such a case and as long as the energy diffusion caused by the collisions is ne- 
glected, which is a reasonable approximation, the wave power can not influence the electron 
relaxation, since the particle never reaches the wave diffusion domain. The relaxation path 
of an electron starting from the LH quasilinear domain (in region 2) can be seen on Fig. 01 
In this case, the supplemental parallel diffusion and drift induced by the LH wave clearly 
lengthen the path, and will be likely to slow the relaxation process, thus enhancing the car- 
ried current. Another possibility is for the initial velocity to belong to region 2, but outside 
the LH domain (m||o < A sample trajectory corresponding to this case is shown on 

Fig. Eland interestingly enough, although the particle has an initial parallel velocity such as 
u\\Q < u\\^i, it experiences the wave influence due to the pitch-angle scattering effect. This is 
the reason why uq appears to be as crucial as M||o. Finally, on Fig. El the relaxation of an 
electron having n||o > 'U||,2 is shown (region 3). Although the wave domain is encountered, 
the net effect on the relaxation length is more complicated, as the upper-velocity boundary 
of the LH quasilinear coefficient, at least in the model chosen here to describe the wave, 
induces a drift towards lower velocities (since ddih/du\\\u^^ ^ < 0) which can accelerate the 
relaxation with respect to a purely collisional trajectory. 

For a quantitative evaluation of the effects of LH waves on the particles thermalization, 
a statistical analysis has to be performed. This is done by considering many electrons with 
initial velocity uq, tracking the relaxation trajectories until thermalization and performing 
the ensemble average, to compute the response function which, in the framework of this 
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stochastic description and according to Eq. UHl is given by 

poo 

Xi(uo)= / dT'{u\\){T') (46) 



When the colhsions are the only effect taken into account, it is possible to average the 
Langevin equations analytically and deduce the response function l|. For our purpose, 
however, this operation is not possible and the computation of the response function has to 
be performed numerically. 

Here, we study {u\\) as a function of time for various initial velocities and values of -Di/i,o, 
the average being performed over 20,000 electrons. The result appears on Fig. |Hl for initial 
positions uq = {u\\q,u±_q) = (3,0), (4,0) and (6,0), and normalized LH diffusion coefficients 
A/i,o = (colhsions only), Ah,o = 0.1, Dih^ = 0.2, Dih^ = 0.4. The colhsional Fisch- 
Boozer solution, given (for Zi = 1) by (n||)(r) = M||o(1 — Sr/nj]) also appears but is perfectly 
superimposed with its numerical counterpart and is thus barely visible. 

It can be observed that for {u\\o,uj_o) = (3,0) or (4,0), the wave clearly delays the 
thermalization. In the third case {up = 6), the effect of the drift induced at the high- 
velocity boundary appears clearly, since the electrons begin by experiencing a faster decrease 
in parallel velocity, on average. However, the energy range of the electrons is largely spread 
by the wave and a significant proportion has not yet thermalized, well after the purely 
collisional relaxation is over. The response function is determined by the balance between 
these two effects and its features shall be presented in the next section. 



V. RESPONSE FUNCTION OF A LHCD PLASMA 

In section IIV[ a numerical method has been presented and employed to perform a ba- 
sic analysis of the velocity space structure when the presence of LH power influences the 
dynamics underlying the current drive process. Although it can be used to compute the 
response function, the perturbation method presented in section ITTTl is more economical in 
terms of computational resources and is thus more adapted to a systematic study of the 
response function properties, which is the goal of the present section. 

In what follows, the same parameters as in section HVl are considered for the Lower Hybrid 
quasilinear domain. On Fig. [7J some level curves the total response function xo + ci-re 
represented in {u\\,u±) space for DihQ = 0.1 (a) and A/i,o = 0.2 (b). For comparison, the 
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corresponding contours of the Fisch-Boozer response function xo appear as dashed hues. 

Several observations can be made about this figure. Firstly, the overall modification of 
the response function is rather moderate, which is consistent with the approximation of 
the adjoint calculation, which requires the collisions to dominate the electron relaxation. 
Secondly, the response function is modified mainly in the LH quasilinear domain, but not 
only and it can be seen to extend to all velocities such as m > u\\^i as well as beyond u\\ = u\\^2, 
which is consistent with the conclusions drawn in section HVl This behavior is clearly visible 
on Fig. IHl where the response function perturbation 6x is shown as a function of /i for various 
values of u, and Dihfi = 0.1. 

This figure shows the strong asymmetrical shape of 6x- As a result, the total response 
function Xi is largely enhanced in the u\\ > region of velocity space, under the influence 
of the LH wave. As predicted for u < u\\^i, we obtain 6x = 0, in other words xi reduces 
to the Fisch-Boozer response function. The pitch-angle scattering effect can cause Sx{u, fi) 
to be non-zero even for u\\ < u\\^i and particularly for u\\ < 0. Another observation is that 
6x falls off rapidly for u > u\\^2 and can even become negative. This effect is best viewed 
when 6x is represented as a function of u\\, for various values of u±, as shown on Fig. IHl For 
completeness, the result from the Langevin equations computation also appears for u± = 0, 
the error bars being deduced from the estimated statistical error. 

One can notice that for m > 6, 5% is negative, which would indicate a deleterious effect 
of the LH wave on the EC-driven current when the latter is carried by electrons excited 
in this region. Although supported by the numerical solution of the Langevin equations, 
this conclusion must, however, be tempered by several considerations: (i) In this region of 
velocity space, Sx is sensitive to the value chosen for A-uy 2 which is not readily available, (ii) 
EC wave absorption at velocities significantly above the upper bound of the LH quasilinear 
domain is difficult and would most likely be impossible in the absence of LH wave. The 
concept of EC current drive improvement or degradation is thus of little sense, (in) The 
Fisch-Boozer response function is proportional to u^, which makes the LH-induced modi- 
fication rather weak for large values of u\\. The latter point is supported by the result of 
Fig. ^[ where Xo and Xo + ^-^e shown as a function of u\\ for the same parameters as 
Fig. El 

This figure shows that beyond the upper boundary of the LH domain, the effect is indeed 
small. Physically, this simply means that no EC efficiency improvement takes place for 
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""11 ~ W||.2 5 aside from the fact that the very presence of a LH plateau is responsible for the 
EC wave absorption at this locationj^. 

More important than the response function itself, as far as the RF current drive is con- 
cerned, is its velocity-space gradient, as is apparent from Eq. E| For EC waves, the differ- 
entiation is to be performed along u±. On Fig. [TTJ the quantity 6x' = d6x/du± is shown 
versus u\\ for various values of u±. 

Whenever u\\ < u\\^2, ^x' appears to be positive, which implies a favorable contribution of 
the LH wave to the EC current. Moreover, even for moderate values of u±, where electrons 
are most easily driven by EC waves, 6x' can be fairly large, provided an appropriate range of 
parallel velocities is selected, which is possible through the use of suitable launching anglesQ]. 
For u\\ ^ U||_25 ^x' can have a negative value, although as stated above, this feature should be 
pondered cautiously. This figure confirms that the LH wave has an overall beneficial effect 
on the current driven EC wave, and that a synergy between the two waves can be expected, 
especally whe„ the latter are excited „r the vein.ty of „,„ < u,,. Note that this parameter 
is generally simply determined by the LH wave accessibility condition [2^, and is therefore 
readily available from experimental measurements of the major plasma parameters. 

VI. ECCD EFFICIENCY IN THE PRESENCE OF LH WAVES 

The features of the distribution function in the presence of LH waves presented in the 
previous section have consequences in terms of ECCD efficiency. From Eq. the total 
current appears as j = jo + ji where jo (Eq. fTI)|) is driven by the LH wave, and is implicitly 
assumed to be unaffected by the presence of the EC wave, in the present model. In section 
mil the response function of the plasma in the presence of LH waves was linearized according 
to Xi = Xo + ^X: which implies that ji, given by Eq.lT^can be cast into the form ji = jec + Sj, 
with 



jec is the EC current obtained when the cross-effect of both waves is not accounted for 
and 6j is a supplemental current, which qualifies for the denomination synergy current (or 
anti-synergy, in the event that 6j and jec have opposite signs). This clear distinction stems 




(47) 



and 




(48) 
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from the linearization introduced in the computation of the response function and allows 
a straightforward separation of the contribution of each process, which is generally the key 
difficulty encountered when trying to characterize a synergy effect. 

Defining jy^ the amount of current generated by wave power p^, the steady-state current 
drive efficiency can be expressed as|jy|] 

dx 



7 / '^"^'^ du 
Vn.^- = -y^ (49) 



au \ 2 

where x is the associated response function and the quasilinear flux. Uf the interaction 
is supposed to be localized in velocity-space, as is generally the case with EC waves and, 
although to a lesser extent, with LH waves also, the efficiency can be approximated to give 



dvi\2 

Evaluating this expression only requires the direction of S^, which unlike its magnitude, is 
only weakly dependent on the precise shape of the distribution function and is well-known^ 
(Sec oc e_L and S^^ oc ey). 

Here, the efficiencies of both waves are computed using Eq. EOwith the same parameters 
as in Section ITTTl A/i,o = 0.1, u\\^i = 3, u\\^2 = 5, A-Uji i = 0.5 and Am||^2 = 1- On Fig. IT^ 
these efficiencies are shown as functions of u\\ for Zi = 1, u± = and u± = 2 

The 4:3 ratio between the LH and EC efficiencies is recovered, in the absence of a cross- 
effectpj]. If this effect is included, the corrected EC efficiency (i.e. (jec + Sj)/Pec) is sig- 
nificantly enhanced in the region of velocity space corresponding to the LH superthermal 
plateau. For the chosen value of Dih^o, it becomes comparable to the LH efficiency. A slight 
anti-synergy effect is observed for u\\ above M||^2 but this point has been extensively addressed 
in Section IIIII The EC efficiency enhancement increases with Uj_ and can exceed the LH 
efficiency on a significant range of parallel velocities. Moreover, for some parameters, the 
counter-current drive amount can be slightly lowered by the synergy effect, which can fur- 
ther enhance the current driven in a EC downshift scheme j^, where the wave interacts with 
electrons having u\\ < as well as My > 0. 

To characterize the ECCD efficiency enhancement, following the definition introduced 
in Ref. 3, the improvement factor Fgyn = {jih+ec ^ jih)/jec is studied. According to the 
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linearization introduced in this model, it simplifies to give 

F„„ = l+«=l + f^ (51) 

3ec dxo/dui_ 

The same parameters as above are used, but the ion charge Zi is varied from 1 to 3. The 
variation of the synergy factor as a function of u\\ is shown on Fig. ^| 

It is seen that the interaction with electrons whose parallel velocity lies in the vicinity of 
the LH quasilinear domain can be very beneficial for the EC current drive. For the chosen 
value of Dihfi, an improvement of the efficiency as high as 40% can be obtained. The anti- 
synergy effect underlined above proves to be marginal. Another observation is that even 
though the plasma ion charge increase is known to be detrimental to superthermal electrons- 
based current drive schemes, it has only a minor influence on the synergy mechanism itself. 
Finally, a particularly noticeable feature is that the improvement factor is found to be weakly 
dependent on the particular velocity space location under consideration. 

Hitherto, the presented simulations have all been performed for Dihfl = 0.1, in order to 
ensure the validity of the perturbation technique employed to derive the response function 
from the adjoint equation. The Langevin equations formalism, introduced in section llVt is 
more demanding from a computational point of view. Still, it has the advantage of offering 
more flexibility than a Fokker-Planck treatment and allows to study how the EC Current 
Drive efficiency depends on Dih^, since it is derived without any assumption regarding its 
particular value. 

This study is performed by solving the Langevin equations (Eq. for 20,000 electrons 
at each initial velocity. For increasing values of Di^ q, the timestep Ar is decreased so 
as to remain small with respect to the variations of Dih (see Eq. l39|) and Fih (Eq. HT| . 
To ensure convergence, however, the computation is performed with several values of Ar, 
and the resulting averaged quantities are extrapolated to Ar 0. For large values of 
Dihfi, this procedure proves to be necessary to avoid the inherent bias induced by the time 
discretization, which adds to the purely statistical errorj27|. 

The same plasma and LH wave parameters as above are chosen, with Zi = 1. The 
simulations are performed for velocity space location {u\\,u±) = (4,1) and the quantity 
under study is 6x = Xi^ Xo where xi is obtained from the Langevin equations (Eq. HHI and 
Xo is given by Eq. 123 As discussed in section IVIf the most relevant quantity, as far as the 
ECCD efficiency is concerned, is 6x' = d6x/du±. To obtain it, simulations are performed 
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for several values of u± and the derivative is obtained numerically. 

On Fig. 1131 Sx and 6x' are plotted versus Dih^, with associated fitting curves. For 
comparison, the adjoint solution for 6x appears also. 

The adjoint solution appears to give a fair result for Di^fl ^ 1, in accordance with the 
validity range of the associated method. Its linear dependence on Di^ q (see Eq.lH^. however, 
leads to an overestimate of the response function, as 6x is found to level off when Dih^o is 
increased. To extrapolate the results as Dihfi tends to infinity, i.e. in a perfectly saturated 
situation, the following fitting function is employed for 6x' 

5x' = ao ■ [1 - eM-aiD^o)) (52) 

A least square fit leads to ~ 7.3, ai ~ 0.7 and 02 ~ 0.5. This gives the extrapolated 
value lim^j^ g^oo i^x' ~ 7.3, or when used in Eq. EH Fgyn ~ 2.8. This means that for the 
parameters considered here, the ECCD efficiency is nearly tripled when compared to its 
"standard" value, i.e. the value obtained in the absence of Lower Hybrid waves. 

VII. CONCLUSIONS 

Owing to the complexity of a full kinetic treatment of the current drive problem in toka- 
maks, added to the difficulty of separating the contributions from various physical processes, 
the existence of a synergy between LH and EC waves has often been disputed. In this pa- 
per, a different approach has been employed to address this question from a theoretical 
standpoint, when the dominant source of deformation of the distribution function is LH 
power. Two complementary methods have been employed: a perturbation solution of the 
adjoint equation allows a fast derivation of the response function including both collisional 
and LH effects in the limit where the former dominate. On the other hand, solving the 
Langevin equations allows to overcome this restriction and in addition to the computation 
of the response function, provides clear insight in the dynamics underlying the process. How- 
ever, they imply a less straightforward mathematical treatment and higher computational 
requirements. 

By application of these two formalisms, it has been shown that a synergy was indeed 
possible between the two waves, provided the EC parameters are chosen to drive electrons 
within or close to the LH quasilinear domain. Even for moderate values of the LH quasilinear 
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diffusion coefficient, a significant improvement of the ECCD efficiency has been obtained. 
Moreover, for sufficiently high values of this coefficient, i.e. when quasilinear saturation is 
reached, the improvement factor appears to be nearly constant. Also, in this study, the 
efficiency enhancement has been found to exhibit only a mild dependence on the particular 
velocity space location. These trends suggest that the synergy mechanism is fairly robust 
and should manifest itself provided the EC waves are launched using a set of parameters 
compatible with the LH quasilinear domain properties. 

In present experiments, the characterization of such an enhancement of EC current drive 
efficiency can be a daunting task. This is mostly due to the fact that, as in numerical 
simulations, a lot of phenomena are involved in the process and are difficult to separate. 
It is nonetheless possible to envision experimental scenarios aimed at studying the LH-EC 
synergy. One such scenario could be creating a fully non-inductive plasma by relying on the 
LHCD system. After a delay equivalent to several resistive times, injecting two EC beams 
with opposite parallel spectra - i.e. opposite toroidal angle, to lowest order - should result 
in zero net EC current in the absence of a cross-effect. The measurement of an additional 
amount of current would then be the signature of a synergy between the two waves. 
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FIG. 1: Velocity-space configuration in the presence of LH waves. Three domains are distinguished, 
each corresponding to a different situation with respect to the wave-induced dynamics. Region (1) 
is such as uo < Region (2) encompasses uq > u^\^i and uyo < U||_2) ^^id Region (3) corresponds 

to U\\o > Ul2- 
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FIG. 2: Sample relaxation path, (a) Collisions only, (b) Collisions + LH wave. The initial position 
{u\'^o,u±o) = (2.9,0), is denoted by an arrow. On (b), the dashed lines delimit the LH domain 
boundaries. 
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FIG. 3: Same as Fig. |^but for initial position (ii||o,^i±o) = (4,0). 
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FIG. 4: Same as Fig. |^but for initial position (ti||o,^*±o) = (2,4). 
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FIG. 5: Same as Fig. |^but for initial position (ii||o,^i±o) = (6,0). 
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FIG. 6: averaged over 20,000 particles as a function of time for u±o = and (a) nyg = 3, 

(b) U||o = 4, and (c) uyo = 5. A/1,0 = (dotted line), Dih^o = 0.1 (dot-dashed line), -D/h,o = 0.2 
(dashed line), Di^ q = 0.4 (solid line). The analytical curve for collisions only appears as a thick 
grayed line. 
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FIG. 7: Iso-contours of the response function xo + Sx a plasma with LH waves such as um = 3, 
liji 2 = 5, = 0.5, Auji 2 = 1 and (a) Di^ q = 0.1, (b) Di^ q = 0.2. The dashed vertical 

lines delimit the LH domain boundaries and the dotted contours represent the collisional response 
function xo- Contours start at |x| = 25 and are equally spaced with |Ax| = 15. 
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FIG. 9: Sx as a function of normalized parallel velocity uy for u± = (solid), 4 (dashed), and 5.5 
(dot-daslied) . Also shown is the response function obtained by numerical solution of the Langevin 
equations for u± = (crosses). 
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FIG. 10: Comparison of the Fisch-Boozer response function xo (thin Hnes) and the response 
function including LH wave effects xo + (thick hnes) for Dihfi = 0.1 for u± = (solid), 4 
(dashed) and 5.5 (dot-dashed). 
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FIG. 12: Current drive efficiency as a function of uy for (a) u± = and (b) u± = 2 for LH waves 
(dashed line) for EC waves alone (thin solid line) and for EC waves when the synergy effect with 
LH waves is taken into account with Dihfi = 0.1 (thick solid line). 
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FIG. 13: Synergy factor in the presence of LH waves for Di/i,o = 0.1 and various values of the 
plasma ion charge: Zi = \ (solid), = 2 (dashed), and = 3 (dot-dashed), (a) u_l = 0, (b) 
U_L = 2. 
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FIG. 14: (a) 8x and (b) d5x/duj_ versus normalized LH quasilinear diffusion coefficient Dihfi for 
U|| = 4 and u± = 1. The solid lines are fitting curves and on (a), the result of the adjoint method 
is shown as a dashed line in the inset containing Sx plotted on logarithm axes. 



34 



